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Abstract. - We use a non-perturbative renormalization-group technique to study interacting 
bosons at zero temperature. Our approach reveals the instability of the Bogoliubov fixed point 
when d < 3 and yields the exact infrared behavior in all dimensions d > 1 within a rather 
simple theoretical framework. It also enables to compute the low-energy properties in terms of 
the parameters of a microscopic model. In one-dimension and for not too strong interactions, it 
yields a good picture of the Luttinger-liquid behavior of the superfiuid phase. 
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Introduction. — In spite of the success of the Bo- 
goliubov theory in providing a microscopic explanation 
of superfluidity [1] , a clear understanding of the infrared 
behavior of interacting boson systems at zero tempera- 
ture has remained a challenging theoretical issue until 
very recently. Besides approximations that do not sat- 
isfy the Goldstone-Hugcnholtz-Pines theorem [2.3], first 
attempts to improve the Bogoliubov theory revealed a sin- 
gular perturbation theory plagued by infrared divergences 
due to the presence of the Bose-Einstcin condensate and 
the Goldstone mode [4,5]. These divergences cancel in 
most physical quantities but lead to a vanishing of the 
anomalous self-energy £ an (<z) in the limit q = (q, ui) — > 
although the linear spectrum and therefore the super- 
fluidity are preserved [6-9]. This observation seriously 
called into question the validity of the Bogoliubov the- 
ory, where the linear spectrum relies on a finite value of 
£ an (<Z — * 0) [10]. The physical origin of the vanishing of 
the anomalous self-energy is the divergence of the longi- 
tudinal correlation function which is driven by the gapless 
(transverse) Goldstone mode - a general phenomenon in 
systems with a continuous broken symmetry [11]. 

The infrared behavior of zero-temperature Bose sys- 
tems is now well understood in the modern language 
of renormalization group (RG) [12-15]. Using a field- 
theoretical renormalization-group approach supplemented 



by the Ward identities associated with the gauge symme- 
try, Castellani et al. were able to establish the exact in- 
frared behavior of a zero-temperature Bose system [14,15]. 
Only for d > 3 does the Bogoliubov theory predict the 
correct infrared behavior, whereas the Bogoliubov fixed 
point is found to be unstable for d < 3 even though the 
low-energy mode remains phonon-like with a linear spec- 
trum. In the approach of Refs. [14,15], the low-energy 
behavior of the correlation functions is expressed exactly 
in terms of thermodynamics quantities such as the density, 
the condensate density or the macroscopic sound velocity. 
Despite its very elegant formulation, this approach how- 
ever does not appear to enable an explicit calculation of 
the correlation functions in terms of the parameters of a 
particular microscopic model. Given the present possibil- 
ities to realize low-dimensional and/or strongly correlated 
Bose systems in ultracold atomic gases [16], it would be of 
great interest to have a theoretical framework allowing for 
quantitative predictions that could be tested against the 
experimental results. 

In this Letter, we use a non-perturbative renormali- 
zation-group (NPRG) technique [17-20] to study inter- 
acting bosons at zero temperature. Not only does our 
approach give the exact asymptotic behavior of the cor- 
relation functions in dimensions d > 1 within a rather 
simple theoretical framework free of infrared divergences, 
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but it also enables to explicitely follow the behavior of the 
system from microscopic to macroscopic scales. In one- 
dimension and for not too strong interactions, it yields a 
good picture of the Luttingcr-liquid behavior of the su- 
pcrfiuid phase. NPRG studies of interacting bosons have 
previously been reported both at finite [21] and zero [22] 
temperature. To a large extent, our results arc comple- 
mentary to those of Ref . [22] . 



Non-perturbative RG approach. 

the following action, 



We consider 



S 



dx 



2m 



(1) 

where tp(x) is a bosonic (complex) field, x — (r, r), J dx = 

J dr J d d r. r <G [0, ft] is an imaginary time, ft — > oo the 
inverse temperature, and \x denotes the chemical potential. 
The interaction is assumed to be local in space and the 
model is regularized by a momentum cutoff |q| < A (with 
A — > oo whenever convenient). We take h = ks = 1 
throughout the Letter. 

The basic quantity of interest in the NPRG is the ef- 
fective action r [</>], which is the generating functional of 
the one-particle irreducible (1PI) vertices. It is obtained 
by a Legendre transform of the free energy lnZ[J] com- 
puted in the presence of an external source term Sj = 
J dx[J* (x)4)(x) +c.c] {4> = (ijj)j) [19,22]. To implement 
the RG procedure, we add to the action an infrared regula- 
tor ASr = / dx ip* (x)R(x — x')tp(x') which suppresses the 
fluctuations with q 2 < k 2 . The functional T[4>] = ^k[<P] 
then becomes k dependent and satisfies the exact flow 
equation 



d t T[4>] = ^Tr [d t R 



(2) 



M+i?)- 1 }, 



(2) 



For a supcrfluid Bose system, the simplest choice reads 

ry>] = r min + J dx{zz x <p(x)d r <\>{x) - v<p*{x)d 2 T <t>{x) 

-Z^{ X )Y-4>{ X ) + ^[n{x)-n f), (3) 

where n(x) = \4>(x)\ 2 is the density. Eq. (3) is obtained 
from an expansion to fourth order about the minimum 
|</>(x)| = y/fiQ, where no denotes the condensate density. 
We use a derivative expansion to order 0(d 2 ) [19,23]. 
For k — > oo, the initial conditions arc Z = Z± = 1, 
V = 0, A = g and no = [i/g, and the effective action 
T[(f>] reproduces the Bogoliubov theory. Although V is not 
present in the original action (1), it is always generated 
by the flow equation (2) [22] and plays a crucial role when 
d < 3. The degeneracy of the minimum |0(a;)| = ^/no 
reflects the gauge invariance (i.e. the U(l) symmetry 
ipW(x) ->■ i) { -*\x)e ±la ) of the action (1). A broken- 
symmetry state can be obtained by picking up a particular 
minimum. It is convenient to write <p — {.4>i + *02)/v / 2 in 
terms of two real fields <fii and 02 and to consider the state 
4> = (y/2no, 0) as an example of broken-symmetry state. 
In Fourier space, the corresponding single-particle vertex 
f( 2 ) = -G" 1 (with Gij = -(ipiipj)) then reads 



Vu 2 + Ze q + 2Xn 

— ZZiLU 



ZZ\LO 

Vur 2 + Ze„ 



(4) 



(e q = q 2 /(2m)). The vanishing of T22(q = 0), which is a 
mere consequence of the U(l) symmetry, naturally imple- 
ments the Hugcnholtz-Pines theorem [3,24] in our formal- 
ism. The combination Ano corresponds to the anomalous 
self-energy E an (g = 0). There are two important quanti- 
ties that can be read off from (4), namely the supcrfluid 
density n s and the Goldstone mode velocity c, 



Zno 



Z/2m 



V+{ZZ 1 ) 2 /(2\n ) 



1/2 



(5) 



where we have introduced the flow parameter t = ln(fc/A). 
In Fourier space, the trace in (2) involves a sum over fre- 
quencies and momenta, as well as a trace over the two 
indices of the complex field <fi. [<p] is the second-order 
functional derivative of T[4>] with respect to <fi. Choosing R 
to diverge for k — > oo, all fluctuations are then suppressed 
and the mean-field theory, where the effective action T[<f>] 
reduces to the microscopic action S[(f>], becomes exact. 
Quantum fluctuations arc gradually taken into account 
by decreasing k and making use of (2). For k = 0, T[(p] 
corresponds to the effective action of the original model 
(1) from which we can deduce all 1PI vertices - and in 
particular the single-particle propagator G = — T^ 2 ) -1 - 
as well as the thermodynamic potential. 

The functional differential equation (2) is too compli- 
cated to be solved exactly. For approximate solutions it is 
sufficient to truncate the most general form of T[<f>] [19,22]. 



For k — > oo, one has n s = hq = ii/g and c = y/ nog/m. 
The supcrfluid density is defined in the usual way from 
the stiffness of the system with respect to a twist of the 
phase of the supcrfluid order parameter <f> = (^/2n,o, 0). 
The expression of the velocity c follows from the equation 
detf W(q) = in the limit q -> 0. 

Inserting (3) into (2), we obtain the flow equations 



d t n = —(d + r) + rii)n 



16s 



A 2 + MA + M 2 - uj 2 
D 2 : 



d t ~X = (d - 2 + 2r] + r/i ) A - 16sA 2 

-5A 3 - 3MA 2 + A\11uj 2 - 6AI 2 } + IMuJ 2 - 4M 3 



D 3 



(0) 
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= -n - 16sA 2 n / <j 



(A + B)(3A- B -4Z 2 u 2 ) 



d t Z 2 = (2 + n + 2 m )Z 2 - 16sA 2 n 



D 2 



2{A + B)(Z 2 B + 1) + AZ 2 cu 2 [Z 2 (3A + 5B) + 2] 
6uj 2 (A + B)[Z 2 {A + B) + \\{2Z 2 B + 1) 



D 4 



d t Q = -(d + 2 + r) 1 )Cl + 8s 



A + M 
D 1 



(7) 



where A = 1 + Z 2 lu 2 , B = A + 2AI, D = AB + uj 2 , M = 
An , v^ 1 = 2 d+1 Tr d / 2 T(d/2), s = (v d /d)[l-r)/(d+2)], and 
/ = fdu)/(2Tr). We have introduced the dimensionlcss 
quantities 



no = ZZ\k no, 
ZZ\Z 2 = e k V, 



A : 

f2 



Z^Z^k^X, 



(8) 



as well as n = —dt In Z and r\\ = —dt In Z\, and chosen the 
regulator i?(q 2 ) = Z(e k — e q )#(efc — e q ) [25]. f2 denotes the 
thermodynamic potential per unit volume in the broken- 
symmetry state (f> = (\/2no, 0). For Z 2 = (i.e. V = 0), 
the integrals over u> can be carried out and we reproduce 
the flow equations derived in Ref. [22]. However, the ap- 
proximation V = — or a mere perturbative treatment 
of V - cannot be used for d < 3 as it predicts the wrong 
exponent (2e instead of e = 3 — d) for the divergence of the 
longitudinal correlation function, the wrong lower critical 
dimension (2 instead of 1), and - for V = - an infinite 
velocity for the Goldstone mode. 

Superfluidity with BEC (d > 1). - When d > 1, 
superfluidity is always accompanied by Bosc-Einstein con- 
densation (BEC): linifc^o no = n^ > 0. For d > 3, the Bo- 
goliubov fixed point is stable; all parameters in the effec- 
tive action F{-o remain finite as k — > 0. V* = limfe~»o V, 
although nonzero, gives only a finite correction to the in- 
frared limit of the vertices. This picture changes dramat- 
ically when d < 3. In this case, both Z\ and A are sup- 
pressed as k — * 0, which explains why the anomalous self- 
energy £ an (<7 = 0) = An vanishes in the infrared limit. 
This suppression is logarithmic for d = 3 and powcrlaw- 
likc in lower dimensions. When d > 1, we can use the fact 
that linifc^o M = limfe_>o Z 2 M = oo to analytically obtain 
the asymptotic behavior for k — -> (table 1). A typical RG 
flow in two dimensions is shown in Figs. 1 and 2. 

The suppression of Z\, together with a finite V* , shows 
that the effective action exhibits a space-time SO(d+l) 
symmetry in the infrared limit [22]. This limit is well 
understood and corresponds to the classical 0(2) model 
in d + 1 dimensions. The symmetry can be made ex- 
plicit by the rescaling f = fcr, f = (Zie^ 1 ^/Z 2 )~ 1 T and 
4>{x) = (ZZiy r Z 2 ~k- d ) 1 / 2 (/)(x) 1 whereby the effective ac- 
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Table 1: Asymptotic behavior for k — > (e = 3 — d). The stared 
quantities indicate nonzero fixed-point values. For d > 1, these 
results are obtained analytically from the flow equations (7). 
For d = 1 one obtains approximate fixed points rather than 
true fixed points (see text). 
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Fig. 1: (Color online) Condensate density no and superfluid 
density n 3 vs. —t for d — 2, no(t = 0) = 0.2 and \(t = 0) = 10. 
Here and in the following figures, we use units where A = 1 
and 2m = 1. 



tion becomes 



J dx\z,Z 2 1/2 



t k (j)* (x)df<t>{x) 
P(x)(d 2 f + V 2 )4>(x) + ^[h(x) - n ] 2 }, (9) 



where n Q = \/Z^n and A' = \j\fZ 2 ~. For a typical fre- 
quency u ~ k, the term linear in df becomes sublead- 
ing with respect to the quadratic one. Equivalently, one 
can observe that for d < 3 the Goldstone mode velocity 
reaches the fixed point value c* = (Z* /2mV*) 1 ^ 2 which is 
independent of Z\. Our numerical results for the scaling 
of n' and A' agree with the known results for the Gold- 
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Fig. 2: (Color online) A, Z\ and V vs. —t for d — 2, no(t = Fig. 3: (Color online) Uq, n* and n* vs. X(t = 0) for d = 2 and 
0) = 0.2 and A(f = 0) = 10. n (t = 0) = 0.2. 



stone regime of the classical 0(2) model in d + 1 dimen- 
sions (table 1). The dimensionless coupling A' vanishes 
for d = 3 and flows to a finite value for d < 3. The rela- 
tion 2rjl +r/2 = —2 (with 772 = —dt In Z2) ensures that the 
Goldstone mode velocity reaches a finite value for k — > 0, 
i.e. that the dynamical exponent takes the value z = 1. 

There are three important relations between 1PI ver- 
tices and thermodynamic quantities that our results 
should fulfill, 



71, = 71 = - 



<9/_t ' 



m(dn/dfi) 



1/2 



\n 



driQ 



(10) 



The equality of the superfluid density n s and the density 
?i is a consequence of Galilean invariance at zero temper- 
ature. Since the chemical potential appears only in the 
initial conditions for the effective action T[<j>], derivatives 
with respect to [i can be numerically calculated by solving 
the flow equations for nearby values of fx. The equality be- 
tween the Goldstone mode velocity c and the macroscopic 
sound velocity c s was proved in Ref. [5]. The last relation 
in (10) is a consequence of gauge invariance [15]. The fig- 
ures 3, 4 and 5 show that the symmetry constraints (10), 
despite a good overall agreement, are not strictly enforced 
in our approach. This can be ascribed to the choice of our 
infrared regulator i?(q 2 ) as well as the Ansatz (3) which 
are both incompatible with Galilean invariance. 

Making use of (10), we can rewrite the propagator 
G(q) = -rW- 1 ^) as 



G 22 (<z) = 
Gu(q) 



2mc 2 no 



1 



n oj 2 + c 2 q? 
mc driQ lo 
n dfi uj 2 + c 2 q 2 ' 
1 



2\n 



(11) 



in the infrared limit. In (11), all quantities except A can be 
evaluated at k = 0. Because of the vanishing of A(fc — > 0), 
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Fig. 4: (Color online) Goldstone mode velocity c* and macro- 
scopic sound velocity c* vs. X(t = 0) for d = 2 and no(t = 0) = 
0.2. 
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Fig. 5: (Color online) Zln* /(A*no) and the condensate "com- 
pressibility" (dno/dfj,)* vs. X(t — 0) for d — 2 and no(t = 0) = 
0.2. 



the longitudinal correlation function Gu diverges loga- 
rithmically in three dimensions and as k~ e below [6,7,9]. 
(see table 1). The dependence on q can be restored by 
\fu> 2 + c 2 q 2 . In Ref. [15], the equa- 



evaluating A at A: 
tions (11) were obtained by imposing the Ward identi- 
ties due to gauge invariance and solving a one-loop RG 
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equation for the sole independent coupling in the limit 
k — ► [26]. By a detailed analysis of the structure of the 
perturbation theory to higher order, it was then argued 
that the equations (11) give the exact asymptotic behav- 
ior. We believe that our RG approach, being intrinsically 
non-perturbativc [19], gives further support to this claim. 

Superfluidity without BEC (d = 1). — In one- 
dimension, as a result of the emerging SO (2) symmetry, 
we find that the long-distance physics is described by the 
classical 0(2) model in d + 1 = 2 dimensions [22]. We 
thus expect the system to be in the "low-temperature" 
phase of the Kosterlitz-Thoulcss phase transition. There 
is no BEC as the condensate density no ~ k v vanishes in 
the thermodynamic limit k — ► 0. However, the superfluid 
density n s = Zuq remains finite. This phase is generally 
described as a Luttinger liquid (LL) characterized by the 
Goldstone mode velocity c* and the LL parameter K [27, 
28]. 

It has been shown that the NPRG gives a good descrip- 
tion of the classical 0(2) model [29,30]. In particular the 
low-temperature phase is characterized by an approximate 
line of fixed points where the beta function becomes very 
small and the running of the renormalized order param- 
eter n' (or, equivalently, the phase stiffness) very slow, 
which implies a very large, although not strictly infinite, 
correlation length £. The anomalous exponent r\ depends 
on the (slowly running) order parameter n' and takes its 
largest value ~ 1/4 when the system crosses over to the 
disordered regime (k ~ 





Fig. 6: (Color online) RG trajectories (n s , A',c) in one dimen- 
sion for various initial conditions no(t = 0) and X(t = 0). The 
points correspond to equal steps in t. 

By solving numerically the flow equations (7) in one di- 
mension, we have obtained very similar results. Figs. 6 and 
7 show the flow trajectories in the space (n s ,c, A') for var- 
ious initial conditions no(t = 0) and X(t = 0). The points 
correspond to equal steps in t so that very dense points 
indicate a very slow running. For a sufficiently small ratio 
X(t = 0)/no(i = 0), we find that trajectories rapidly hit 
an approximate plane of fixed points defined by A' ~ 15, 
where the running of the superfluid density n s and the 
Goldstone mode velocity c become very slow. As for the 
classical 0(2) model, we infer from this observation that 




A' 



Fig. 7: (Color online) Same as in Fig. 6 but in the plane (n a , A') 




Fig. 8: (Color online) Anomalous exponent r\ vs. —t for various 
initial conditions no (t = 0) and A(f = 0). 



the correlation length £ is extremely large for these tra- 
jectories. For very long RG time —t (k ~ C" 1 ), the sys- 
tem eventually crosses over to the disordered regime. On 
the approximate plane of fixed points, the scale-dependent 
anomalous exponent r\ varies slowly about a value that de- 
pends both on n s and c (Fig. 8). It then reaches its max- 
imum value 7] c ~ 0.29 - to be compared with the exact 
exponent r} c = 1/4 at the Kosterlitz-Thoulcss transition of 
the classical 0(2) model - before rapidly dropping to zero 
once k ~ For M 3> 1, which corresponds to a large 
superfluid density n s , we obtain the analytic expression 
77 = toc/(27T7i s ) from (7). This in turn determines the LL 
parameter K = 1/ (277) = irn s / (mc), which is the expected 
value in a Galilean invariant system where n s = n [27]. 

Thus for a sufficiently small ratio X(t — 0)/n s (or c/n s ), 
we obtain a good picture of the Luttinger liquid behavior 
of the superfluid phase. When this ratio is too large, the 
flow trajectory does not reach the approximate plane of 
fixed points, and the system is in the "high-temperature" 
(disordered) phase of the classical 0(2) model. This result 
is in contradiction with known results in one dimension 
where the action (1) corresponds to the exactly soluble 
Licb-Liniger model [31,32]. This model is parameterized 
by the dimcnsionless parameter 7 = mX(t = 0)/n and its 
low-energy description is a Luttinger liquid with a parame- 
ter K = if (7) varying in the interval [1, 00 [ as 7 decreases 
from infinity to zero [28]. The limit 7 —> 00 (K = 1) 
corresponds to hard-core bosons. Thus the anomalous ex- 
ponent 77 should take its highest value l/(2K min ) = 1/2 
for 7 — >■ 00 rather than 1/4 as predicted by our results. 
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A possible explanation for the failure of our approach to 
correctly describe the strong-coupling limit of the action 
(1) in one dimension is the derivative expansion used in 
the Ansatz (3) for F[</>]. Quite generally, the derivative 
expansion is known to work best when 77 is small [19,20]. 

Conclusion. — The NPRG technique discussed in 
this Letter provides an efficient method to control the 
infrared divergences appearing in the perturbation the- 
ory of zero-temperature Bose systems. It extends the ap- 
proach of Rcf. [22] and reproduce the results obtained ear- 
lier by a ficld-thcorctical RG approach combined with the 
implementation of Ward identities due to gauge invari- 
ance [14,15]. The non-trivial infrared behavior in dimen- 
sions 1 < d < 3, characterized by the divergence of the 
longitudinal correlation function and the vanishing of the 
anomalous self-energy E an (q — > 0), turns out to be related 
to the emergence of a space-time SO(d+l) symmetry at 
low energy. This implies a close link between the super- 
fluid phase and the Goldstone regime of the classical 0(2) 
model in d + 1 dimension [22] . 

Our approach also describes one-dimensional systems 
where superfluidity exists without BEC in the thermody- 
namic limit. The supcrfluid phase exhibits a Luttinger- 
liquid behavior that is well captured by the NPRG ap- 
proach for weak interactions. Although our results, based 
on a derivative expansion of the effective action T[0] , break 
down at strong coupling, they might be improved by a 
more refined treatment of the momentum dependence of 
the vertices [33]. 

An important feature of the NPRG is that it not only 
yields the infrared behavior of correlation functions but 
can also compute propagators in terms of the parameters 
of a microscopic model. It thus provides an efficient tool 
for the explicit calculation of physical quantities beyond 
the Bogoliubov theory while satisfying basic requirements 
such as the Hugenholtz-Pines theorem as well as yielding 
the correct infrared behavior, a task that has been known 
to be difficult in interacting boson systems [2,34]. 
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